Phase diagram of hard-core bosons on a frustrated zig-zag ladder 



Davide Rossini 

NEST, Scuola Normale Superiore & Istituto di Nanoscienze - CNR, Piazza dei Cavalieri 7, T56126 Pisa, Italy 

Valeria Lante and Alberto Parola 
Dipartimento di Fisica e Matematica, Umversitd dell'Insubria, Via Valleggio 11, 1-22100 Como, Italy 

Federico Becca 

Democritos Simulation Center CNR-IOM Istituto Ojficina dei Materiali, Via Bonomea 265, 1-34136, Trieste, Italy 

(Dated: January 18, 2011) 

We study hard-core bosons with unfrustrated nearest-neighbor hopping t and repulsive interaction 
V on a. zig-zag ladder. As a function of the boson density p and V/t, the ground state displays 
different quantum phases. A standard one-component Tomonaga-Luttinger liquid is stable for p < 
1/3 (and p > 2/3) at any value of V/t. At commensurate densities p — 1/3, 1/2, and 2/3 insulating 
(crystalline) phases are stabilized for a sufficiently large interaction V. For intermediate densities 
1/3 < p < 2/3 and large V/t, the ground state shows a clear evidence of a bound state of two 
bosons, implying gapped single-particle excitations but gapless excitations of boson pairs. Here, 
the low-energy properties may be described by a two-component Tomonaga-Luttinger liquid with 
a finite gap in the antisymmetric sector. Finally, for the same range of boson densities and weak 
interactions, the system is again a one-component Tomonaga-Luttinger liquid with no evidence 
of any breaking of discrete symmetries, in contrast to the frustrated case, where a Z2 symmetry 
breaking has been predicted. 
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I. INTRODUCTION 

Quantum many-body systems can give rise to remark- 
able collective states of matter that have no counterpart 
in their classical analogs. Archetypal examples include 
superfluids, superconductors, and insulating quantum 
liquids that play a major role in modern Condensed Mat- 
ter and Materials Science. However, collective quantum 
phenomena are also ubiquitous in other contexts, like Nu- 
clear Physics, Quantum Chemistry, and Atomic Physics. 
In particular, ultra-cold atoms loaded into optical lattices 
provide a unique possibility for engineering quantum sys- 
tems with a very high degree of tunability and control of 
the experimental parameters.— They allow the realiza- 
tion of "quantum simulators" for ideal Condensed Mat- 
ter models, which may provide answers to fundamental 
questions.^ The first striking demonstration in this direc- 
tion has been the observation of the superfluid to Mott 
insulator transition for bosons with short-range interac- 
tionsr more recently, a fermionic Mott insulator has been 
also observedi^ Theoretical progresses and experimen- 
tal achievements steadily open new research directions. 
Highly non-trivial phenomena and very rich phase dia- 
grams are now conceivable by considering further ingre- 
dients, like long-range interactions, spin or multi-species 
models, frustration, and disorder. 

Cold gases of bosonic particles trapped in optical lat- 
tices may be very well described by simple Bose-Hubbard 
models;^ which contain hopping and short-range interac- 
tion termsi^ii Models of strongly-interacting bosons in 
one-dimensional (ID) or quasi-lD lattices constitute im- 
portant examples where unconventional phases can be 



stabilized at low temperature. These systems do not rep- 
resent a purely abstract problem, since it is now possi- 
ble to confine the atomic species in almost decoupled ID 
tubes, with^ or without^ an optical lattice. In realistic 
experimental setups, beyond pure ID systems, the case 
of a two-leg ladder is very easy to obtain. Indeed, one can 
realize a double-well potential along a direction (say, y) 
like in Ref. 10, and a potential creating a cigar geometry 
in the x-axis. Then, by superimposing a further peri- 
odic potential along x, one realizes a two- leg Hubbard 
model. By playing with the distance between tubes and 
the height of the barrier between the two legs, one could 
tune the hopping rate between the legs. Likewise, the 
intra-chain hopping rate can be tuned by appropriately 
setting the strength of the periodic potential along the x 
direction. 

Recent experimental results have driven a new impetus 
to understand their relevant low-energy properties. In 
particular, the study of quasi-lD systems, e.g., ladders, 
may be very important, in order to elucidate the na- 
ture of exotic quantum phases that escape the standard 
Tomonaga-Luttinger theoryj^i From a theoretical point 
of view, ladders are quasi-one-dimensional systems. For 
this kind of anisotropic systems, very efficient and accu- 
rate numerical methods have been also developed in the 
last 20 years, like exact diagonalizations by the Lanc- 
zos techniquei^, or the density-matrix renormalization 
group (DMRG) method<^ Within these approaches, it is 
possible to have numerically exact results on fairly large 
clusters, so to have insights into the physical properties 
at the thermodynamic limit. Moreover, in ID systems a 
clear theoretical framework is provided by the bosoniza- 
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tion methodr^ which is helpful for classifying the possible 
phases. 

In the following, we will consider a zig-zag ladder with 
two legs, which is topologically equivalent to a ID lat- 
tice with equal nearest and next-nearest neighbor hop- 
ping and interaction, see Fig. [TJ We will study the 
case of hard-core bosons that interact through a nearest- 
neighbor potential: 

'H = -tY,{blbj+h.c.) + V^n,nj, (1) 

where (i, j) indicate nearest-neighbor sites in the zig-zag 
geometry of Fig. [3 b] (bi) creates (destroys) a boson on 
the site i and rii — h\hi is the boson density. The hard- 
core constraint is imposed via the additional requirement 
Ui — 0,1 on each site. The number of sites and bosons 
will be denoted by L and M, respectively. The boson 
density will be denoted hj p = M/L. In the following, 
we will focus on the case of unfrustrated hopping, i.e., 
t > 0. 

We mention that the Hamiltonian ([1]) maps onto a sys- 
tem of 5* = 1/2 spins {S^ — rii ~ 1/2) with antiferro- 
magnetic coupling — V between the z components 
of spins and super-exchange J^^ — ~2t between their x 
and y components. Previous works have been focused 
on the case with negative hoppings, i.e., t < Q, that cor- 
responds to an antiferromagnetic Heisenberg model. In 
particular, the one-dimensional chain with SU(2) sym- 
metry (i.e., V — 2\t\) but different nearest- (Ji) and 
next-nearest-neighbor ( J2) interactions has been widely 
discussed, also in presence of a finite magnetic field J^"— 
On the other hand, here we are interested in the case with 
positive hopping parameters and V > 2t, in order to de- 
scribe strongly interacting bosons in low-dimensional sys- 
tems that are relevant for atomic gases trapped in optical 
lattices. However, we will show that some features of the 
phase diagram do not depend upon the sign of t and can 
be understood on the basis of the strong-coupling (clas- 
sical) limit V/t — 7> 00. 

For future reference, it is useful to point out that in lat- 
tices where each site has the same coordination number 
z (equal to the number of neighbors), the microscopic 
model is invariant under "particle-hole symmetry" , de- 
fined as the spin rotation of an angle tt around the x- 
axis. This particle-hole transformation, which is indeed a 
canonical transformation due to the hard-core constraint, 
allows one to relate the full energy spectrum at different 
densities: 

Vz 

E{L - M) = E{M) + — {L- 2M) (2) 

giving rise to a phase diagram symmetrical across the line 
p — 1/2. Open-boundary conditions violate this symme- 
try leading to specific finite-size effects, since a few sites 
have a smaller coordination number. In the following, 
open-boundary conditions will be used in DMRG, while 
periodic boundary conditions will be adopted in Lanczos 
calculations. 
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FIG. 1: The two-leg ladder (a) is topologically equivalent to 
a one-dimensional chain with nearest and next-nearest con- 
nections (b). 
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FIG. 2: (Color online) Examples of classical ground states 
for different boson densities. Full (empty) circles indicate 
particles (empty sites), (a) p — 1/3, where the ground state is 
three- fold degenerated; (b) 1/3 < p < 1/2, where the ground 
state has a finite entropy; (c) p = 1/2, where the ground state 
is four-fold degenerated and corresponds to pairs of nearest- 
neighbor particles. 

The paper is organized as follows: in Sec. |TT1 we il- 
lustrate some useful preliminary consideration related to 
the classical limit of i = 0; in Sec. IIIIl we present our 
numerical results; finally, in Sec. lIVl we draw our conclu- 
sions. 



II. PRELIMINARY CONSIDERATIONS 

Before presenting our numerical results, we would like 
to discuss in some detail the classical limit oi t = 0, 
which is relevant for the physical properties in the strong- 
couphng regime, i.e., V/t 1. Let us take for conve- 
nience L = 3 X n with n multiple of 4 and quantize the 
q vectors as g = ^k, with fc = 0, . . . , i — 1. Whenever 
the number of bosons M < n, i.e., p — M/L < 1/3, 
the classical ground state has zero energy and finite 
entropy. A finite value of t/V will drive the system 
towards a Tomonaga-Luttinger liquid with long-range 
charge-density-waves or superfluid correlations. Exactly 
at M = n (namely p = 1/3) a commensurate density 
wave of longitudinal wave-vector q = 2n/3 sets in, lead- 
ing to a three-fold degenerate solid-like ground state, as 
shown in Fig. [2Ja). The energy gap is proportional to V 
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and, therefore, this insulating state is expected to remain 
stable also in presence of a small hopping parameter (see 
below, our numerical results). By further adding bosons, 
we find that for any density of the form M — n + 2S, 
with S equal to an integer smaller than rt/4 (so that 
1/3 < p < 1/2), the ground state is highly degenerate 
and still gapped, with energy E = SVS and single par- 
ticle gap E{M + 1) - E{M) = 2V. All possible ground 
states may be obtained by viewing the state, in the one 
dimensional topology, as a mixture of single particles and 
nearest-neighbor pairs, and placing these objects on the 
lattice in such a way that the nearest-neighbor sites of 
each object are empty, see Fig. HJb). Given the huge 
degeneracy of all these classical states, a finite value of 
t has a dramatic effect, as it will be shown in the next 
section. Note that for these densities the binding energy 
A = E{M) + E{M + 2) -2E{M + 1) is negative (i.e., 
A — ~V) implying the formation of boson pairs in the 
model. Finally, for p = 1/2, only nearest-neighbor pairs 
are present and the ground state is a four-fold degen- 
erate gapped "molecular solid" characterized by a crys- 
tal ordering with longitudinal wave-vector q = 7r/2, see 
Fig. He). 

These results for the two-leg ladder in the classical 
limit can be conveniently summarized by considering the 
boson density p as a function of the chemical potential p,: 
at /i = 0, the density jumps from p = to p = 1/3, where 
it displays a plateau up to p — 3X^/2. Then, the density 
jumps again to p = 1/2 and then remains constant up 
to p = bV/2. A density discontinuity from p — 1/2 to 
p — 2/3 is followed by a plateau to p — AV where p 
jumps again to its limiting value p = 1. Therefore, in the 
classical limit, the model displays three distinct gapped 
(sohd) phases of density p ~ 1/3, 1/2, and 2/3, besides 
the trivial "empty" (p = 0) and "full" (p = 1) states. At 
intermediate densities phase coexistence between neigh- 
boring phases sets in. As we will see in the next para- 
graph, a finite value of t removes this coexistence by fa- 
voring stable phases. In particular, for 1/3 < p < 1/2, 
long-range charge-density-wave correlations develop at 
an incommensurate wave-vector g, which smoothly in- 
terpolates between the two limiting values (g = 27r/3 for 
p = 1/3, and g = 7r/2 for p = 1/2) characterizing the two 
solids at coexistence for t — 0. 
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FIG. 3: (Color online) Boson density p as a function of the 
chemical potential p for three different values of the inter- 
action V/t. Calculations have been done on a lattice with 
L — 108 sites and open boundary conditions. 
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FIG. 4: (Color online) Entanglement entropy S as a function 
of the (reduced) block length x, for p = 5/12 and for differ- 
ent values of the interaction strength V/t. The dashed line 
indicates the slope 1/6 (equivalent to c = 1). The number of 
sites is L = 216. 



III. RESULTS 

Here, we present our numerical results on the two- 
leg ladder. The standard finite-size DMRG algorithm 
has been adopted)^ fixing the total number of sites L 
and bosons M. Systems of up to L = 216 sites have 
been simulated, keeping m ~ 300 states and performing 
nsw 6 sweeps, in order to reach convergence both in 
the energies and in the measure of observables. Open- 
boundary conditions have been used, see Fig. [IJb). The 
latter choice determines small asymmetries in the phase 
diagram when comparing p and (1 — p), which vanish in 



the thermodynamic limit. In addition, DMRG results 
have been compared to Lanczos diagonalizations up to 
L = 36 sites with periodic boundary conditions. 

In Fig. [3l we report the behavior of the boson den- 
sity p as a function of the chemical potential p: for each 
p, the corresponding density p = M/L is obtained by 
minimizing the free energy E{M) — pM with respect 
to M. For small interactions, on any finite size, p{p) 
displays small steps for every value of M; therefore, 
the ground state is always gapless, the density being a 
smooth function of the chemical potential p. For p < 1/3 
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FIG. 5: (Color online) Difference between the energy needed 
to add a particle /i^ = E{M + 1) — E{M) and the one needed 
to remove it n~ = E{M) - E{M - 1) (where E{M) is the 
energy of M particles) as a function of the density p, for 
different sizes of a system with V/t = 5 (upper panel) and 
V/t — 10 (lower panel). 




FIG. 6: (Color online) Boson density p as a function of the 
chemical potential fi for positive and negative values of the 
hopping parameters. Calculations have been done on a sys- 
tem with L = 108 sites with open boundary conditions. The 
existence of a plateau at p = 1/2 is clear for t — —1, whereas 
no plateau is present for t — 1. 



(and p > 2/3), the system can be described by a one- 
component Tomonaga-Luttinger liquid, where the low- 
energy excitations are free massless bosons ((/), 0) and its 
central charge is c = iJ^ The latter quantity can be nu- 
merically measured through the entanglement entropy: 



S{1) = -Tro [prdmil] InprdmCO] 



(3) 



where PrdmiJ) is the reduced density matrix for the sub- 
system ri, containing all the sites from 1 to I, and is 



p,d,^(Z) = Trj-,|*o)(*o|. 



(4) 



Here |\l/o) is the ground state wave function and defines 
the environment (all sites from / -|- 1 to L). The central 
charge can be obtained from the logarithmic divergence 



6 



(5) 



X — In 


■ 1 




— sm 






TT ' 





which is valid in the thermodynamic limit and / ^ 1. In 
finite systems, it is useful to consider)^ 



(6) 



instead of InL Since the reduced density matrix Prdm{l) 
is directly available by DMRG calculations, the central 
charge c can be easily computed, by fitting the linear 
slope of S{1) as a function of x. We verified that c = 1 
in the low-density regime p < 1/3 and at any interaction 
strength. 

The intermediate boson density at weak coupling re- 
quires a deeper discussion. Indeed, a naive analysis based 
upon the band structure (and the connection to spin-less 
fermions through the Jordan- Wigner transformatiorti^) 
would suggest that the low-energy physics is described 
by a two-component Tomonaga-Luttinger liquid, with 
c = 2. However, we find that c = 1 also in this regime, 
see Fig. HI The limiting case with V — Q corresponds 
to "non-interacting" hard-core bosons with unfrustrated 
hopping amplitude, which are expected to give rise to 
a standard quasi-condensed quantum liquid. Therefore, 
spin-less fermions and hard-core bosons are inherently 
different in this regime of densities (e.g., the Jordan- 
Wigner string caused by hopping along the legs plays 
a relevant role in this zig-zag geometry). We mention 
that in the case with frustrated hopping (equivalent to 
the antiferromagnetic spin XY model), the ground state 
is expected to develop a chiral order that breaks a dis- 
crete Z2 symmetry.— 1^ Similarly, a chiral phase has been 
also predicted for isotropic frustrated spin chains in pres- 
ence of a magnetic field. ^^'^^'■^'^ On the contrary, for our 
unfrustrated model, we do not find any evidence of sym- 
metry breaking for y > 0, as can be seen from Lanc- 
zos spectra at low-energy (not shown). Therefore, we 
conclude that there is a fundamental difference between 
frustrated and unfrustrated bosons at these densities: for 
the former case, a gapless state with broken Z2 symme- 
try is expected, while, in the latter one, a pure gapless 
Tomonaga-Luttinger liquid is realized. 

By increasing the ratio V/t, two plateaus emerge at 
p = 1/3 and 2/3, for V/t > 8, indicating the stabiliza- 
tion of insulating phases with a finite excitation gap, see 
Fig. [21 These states can be adiabatically connected to 
the solid phases that have been discussed in Sec. HT] Re- 
markably, the stabilization of crystalline phases at com- 
mensurate fillings is accompanied by a significant mod- 
ification of the intermediate phase at 1/3 < p < 2/3. 



FIG. 7: (Color online) Density structure factor N{q) for differ- 
ent boson densities. Calculations have been done on lattices 
with L — i08 (solid curves) and L = 60 (dashed curves) for 
V/t = 20. 



Indeed, whereas for V/t < 8 the liquid phase displays 
a standard gapless behavior, with a vanishing excita- 
tion energy when adding or removing a single boson, 
the strong-coupling phase is instead gapped for single- 
particle excitations, as clearly shown in Fig. [Sj where 
the difference between ^+ = E{M + 1) - E{M) and 
^- = E{M) - E{M - 1) is reported. In the strong- 
coupling regime, for p < 1/3, /i^ — /i^ is a smooth 
function of p that goes to zero in the thermodynamic 
limit, while for intermediate densities, this quantity re- 
mains finite, indicating the presence of a finite gap in 
the single-particle excitations. However, such a strong- 
coupling phase is not insulating, since excitations of pairs 
of bosons are still gapless. The plot of p{p) shown in 
Fig. [3] supports this interpretation showing that, on any 
finite-size system, steps twice as big as in the standard 
Tomonaga-Luttinger liquid occur. The emergence of the 
two-boson bound state is simply due to a potential energy 
gain and does not depend on the sign of t: as previously 
noticed, even in the classical limit {t — 0), the binding 
energy A = E{M) + E{M + 2)-2E{M+l) is negative in 
this part of the phase diagram (see Sec.|ll|. Not surpris- 
ingly, a similar behavior has been also found in the J1—J2 
Heisenberg model in presence of a magnetic field] ^^'^^1^^ 
In spin systems, the effect, which was called "even-odd" , 
is characterized by the existence of two-magnon excita- 
tions in the region of weakly coupled antiferromagnetic 
chains. This phase may be described by a low-energy 
theory with two bosonic fields {(f>m9n), with n = 1,2, 
which give rise to symmetric and antisymmetric combi- 
nations; the presence of relevant interactions leads to the 
opening of an energy gap in the antisymmetric channel, 
implying a two-boson bound stateJ^ii^ 

We would like to mention that the solid at p = 1/2 ap- 
pears to be less stable than the other two at 1/3 and 2/3. 
Indeed, we start to detect its existence for V/t = 10-^ 11, 
although the width of the plateau in a p{p) plot is still 



FIG. 8: (Color online) Position of the peak q of the density- 
density correlations as a function of the density. The value of 
q does not depend upon the interaction strength V/t. 




FIG. 9: (Color online) Tomonaga-Luttinger parameter K as 
a function of the boson density p for V/t = 10. 



fairly small. In Fig. [Sj we show that eventually a clearly 
insulating phase exists for V/t — 20. On the contrary, a 
very stable solid phase for p = 1/2 may be obtained by 
changing sign of the hopping parameters; in this case, the 
model is equivalent to the anisotropic J1—J2 Heisenberg 
chain, which shows a dimerized phase at zero magnetiza- 
tioni^iS£ In Fig. [51 we report a comparison of the curves 
p{p) for positive and negative hopping parameters and 
V/\t\ = 1. Only when t — —1 a clear plateau is present, 
indicating that change in the sign of the hopping pa- 
rameters has a dramatic effect in the stabilization of the 
insulating phase at half filling. 

The nature of the ground state can be better character- 
ized by studying the density-density correlation function. 
In the case with open-boundary conditions, used within 
DMRG calculations, we define: 

N{r,r') ^ {nrfir') - {nr){nr'), (7) 

which depends separately on r and r' . It is convenient to 
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FIG. 10: (Color online) DMRG results for the density-density 
correlation functions on systems with open boundary condi- 
tions. Left panel: small-g behavior of N[q) for a 216-site 
system. Data for 36 sites with periodic boundary conditions 
obtained by Lanczos diagonalizations are also shown. The 
boson density is p = 5/12, while V/t = 10. A linear fit of nu- 
merical DMRG data gives K = 0.27 ±0.01, while for Lanczos 
data K = 0.29 ± 0.02. Apart from the system sizes, the small 
discrepancies could be ascribed to different boundary condi- 
tions imposed in the two cases. Right panel: size scaling of 
the peak of the structure factor N{q) oc L" . The exponent 
a = 0.71 ± 0.01 is in good agreement with the value 1 — K 
obtained from the previous fits. Dashed lines are fits of the 
numerical results. 



perform the Fourier transform 
2 



N{q) 



L 



-^sm{qr)sm{qr')N{r,r'), (8) 



with q — lunj (L + 1) and n — 1, . . . , L. With periodic- 
boundary conditions, implemented with Lanczos algo- 
rithm, correlation functions only depend upon the dis- 



tance \r - 



and the standard Fourier transform can be 



employed. In Fig. [3 we show the behavior of the density 
structure factor N{<q) at V jt = 20 for different densities. 
For p < 1/3, the ground state is described by a standard 
one-component Tomonaga-Luttinger liquid, for which (in 
a translationally invariant system): 



K 



A 



cos(27rpr) 

„2K 



(9) 



where A is a numerical constant and K is the Tomonaga- 
Luttinger parameter that describes the low-energy 
modeliii N{q) displays a peak a.t q — 2np (see Fig. [8]), 
in agreement with a single-band picture with kp — irp. 
No long-range order is implied, since the finite-size peak 
grows slower than L: N{q) - L^^'^'^ . The one particle 
density matrix decays as l/r^^^^. Therefore, the dom- 
inant correlations have a charge-density wave (CDW) 
character for K < 1/2, while for if > 1/2 a diver- 
gence in the q — occupation number occurs (i.e., quasi- 
condensation) . 

A convenient way to compute the Tomonaga-Luttinger 
parameter K is via the small-g behavior of N(q). From 
Eq. ^ we easily obtain N{q) = j- q. In Fig. [HI we report 



the resulting if as a function of the density p for V/t = 10 
(also obtained by Lanczos diagonalizations on L = 36 
sites). At low density, bosons are in a quasi-condensed 
state but, close to quarter filling, the crossover to a CDW 
phase takes place. 

At the commensurate density p = 1/3 and sufficiently 
large V/t, true long-range order sets in and N{q) ~ L, 
together with a quadratic behavior of N{q) at small mo- 
menta, i.e., N{q) ^ q^ . When approaching p = 1/3 by 
varying the boson density, the parameter K is expected 
to remain finite, with a value that does not depend upon 
the interaction strength and is twice the limiting value 
obtained at constant density (p = 1/3) when V/t tends 
to the transition value from below We mention that, 
in the classical limit t = Qhy using the results of Ref. H^, 
we find that K = l/9forp— >■ 1/3, which may be com- 
patible with the strong renormalization that is observed 
in our numerical results for V/t — 10 when approaching 
the insulating phase, see Fig. [21 

For intermediate boson densities, i.e., 1/3 < p < 1/2, 
and strong interactions the model is described by a two- 
component boson field with a gap in the antisymmetric 
channel. In this case 



, y 2 K , 
{nrTlo) = p — + A 



^,cos [7r(l -p) r] 



(10) 



where A' is a numerical constant and K is the Tomonaga- 
Luttinger parameter in the effective low-energy model 
of the gapless bosonic field. Pair condensation can be 
investigated trough the correlation function 



B 

A/K 



(11) 



where i? is a numerical constant. In this regime, N{q) = 
^q and the peak of the structure factor shifts to g = 
7r(l — p), see Fig. |51 Again there is no true long-range 
order, since N{q) diverges as L^^^ . By fitting the small- 
q part of the structure factor, we are able to extract the 
behavior of K also in this part of the phase diagram, see 
Fig. [9] for the case of V/t = 10. A fairly good agreement 
is obtained by comparing the power law divergence of 
the peak in N{q) with the expected exponent 1 — K. 
In Fig. [TOl we report these results for p — 5/12 and 
V/t = 10. In this phase, CDW correlations dominate 
and coexist with power-law pairing correlations, which 
however do not lead to (quasi) pair condensation because 
if < 1 in the whole density interval. 

Finally, for p = 1/2 and large V/t, the ground state has 
again crystalline order with N{q) ^ L and N{q) ^ q^. 



IV. CONCLUSIONS 

In this paper, we have studied hard-core bosons with 
unfrustrated hopping and nearest-neighbor interaction 
on a frustrated zig-zag ladder. For small V/t, the ground 
state is a gapless Tomonaga-Luttinger liquid for all den- 
sities. For p < 1/3 (and p > 2/3), the density-density 
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correlations have a peak for q = 2iTp. This fact sug- 
gests that a fermionic band picture (originating from the 
standard Jordan- Wigner transformationi^) can be cor- 
rect and the Jordan- Wigner string is a smah perturba- 
tion to free spin-less fermions. For intermediate densities 
1/3 < p < 2/3, N{q) shows a peak at q — tt{1 — p), 
still having one gapless mode (e.g., c — 1). This fact 
contrasts the naive expectation based upon the Jordan- 
Wigner transformation. Moreover, we do not find any ev- 
idence of discrete Z2 symmetry breaking, as found in the 
case with frustrated hopping.— 1^ In this case, for weak 
interactions, the ground state is a pure quantum liquid, 
with a single gapless mode. By increasing V/t, three solid 
phases appear at commensurate boson densities, first at 
p = 1/3 and 2/3 (for V/t ~ 8) and then at p = 1/2 (for 
V/t c^i 10): these are insulating phases with long-range 
order in the density profile. For 1/3 < p < 2/3 there is 
CDW phase, which has gapped single-particle excitations 
but gapless excitations for pairs of bosons. This phase, 
which has been also found in frustrated antifcrromag- 



netic spin models, can be described by a two-component 
Tomonaga-Luttinger model with a finite gap in the anti- 
symmetric sector. 
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